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Sufficient conditions for the wave instability in general three-component reaction-diffusion systems 
are derived. These conditions are expressed in terms of the Jacobian matrix of the uniform steady 
state of the system, and enable us to determine whether the wave instability can be observed as the 
mobility of one of the species is gradually increased. It is found that the instability can also occur 
if one of the three species does not diffuse. Our results provide a useful criterion for searching wave 
instabilities in reaction-diffusion systems of various origins. 

"cW : 

] The wave instability provides an important mechanism for pattern formation in nonequilibrium chemical systems. 
. When it takes place, a critical mode corresponding to a traveling wave with a certain wavenumber and oscillation 
frequency begins to grow, destabilizing the uniform steady state. Although beiiig less known, the wave instability 
, has already been considered in 1952 by A. Turing in his pioneering publication where the classical (i.e. static) 
■ Turing instability, leading to the establishment of a periodic stationary pattern has also been introduced. Therefore, 
' it may be also appropriate to describe it as the oscillatory Turing bifurcation. Moreover, it was already noticed by 
^ . A. Turing [ij that at least three interacting species arc needed for this instability to occur. 

Because of the spatial reflection symmetry, waves traveling in the left and right directions have the same growth 
rates and both of them begin to spontaneously develop above the instability threshold. Nonlinear interactions between 
^0 ' such modes determine whether one of the modes gets suppressed, so that a wave traveling in a certain direction is 
Ph ■ established, or standing waves, representing superpositions of left and right traveling waves, are instead formed 
(-H ' The wave patterns resulting from such instability can also exhibit further instabilities, and wave turbulence may set 
on. 

In contrast to the classical Turing bifurcation, which has been extensively discussed for both biological and chemical 
systems Jl|, Q-Q; the wave bifurcation has so far attracted less attention. It has been considered for special chemical 
models [1] and its existence was suggested in the experiments with Belousov-Zhabotinsky microemulsions [l^. There 
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' are also publications in which this instability was discussed for special ecological models 
. Because at least three species are needed for the wave instability to occur, the linear stability analysis is more 
00 ' complex in this case, as compared with the classical Turing bifurcation in two-component activator-inhibitor systems. 
, The complexity of the stability analysis, which has been previously performed separately for individual chemical 
' systems, has probably also been responsible for the fact that the wave instability has not been broadly investigated 
' for reaction-diffusion media. 

I In this article, we present a general derivation of the sufficient conditions for the wave (or oscillatory Turing) 
CO ' bifurcation in arbitrary three-component reaction-diffusion models. The final conditions are formulated in terms of 
the elements of the Jacobian matrix of the uniform steady state. They tell whether the wave bifurcation is possible 
when the mobility of any chosen species is gradually increased, while diffusion coefficients of other species are kept 
constant. As we show, the instability is possible even if one of the three species is immobile. Using our general 
expressions, the analysis of the wave bifurcation in specific reaction models will be largely simplified. 
I We consider reaction-diffusion systems with three chemical reactants U, V and W. Local densities of the reactants 
" " are denoted a,s u = [U],v = [V] and w ~ [W] . All reactants diffuse over the space and undergo local chemical 
reactions. Generally, such systems are described by equations 

' du 

— = f{u, V, w) + D„V^u, 

^ = giu,v,w)+DyV^v, (1) 
at 

— = h{u, V, w) + U^V w, 
I dt 

where functions /, g and h represent the local reactions. Diffusion coefficients of the reactants are Du,Dy and D^. 
We assume that a uniform steady state {u,v,w) = {u,v,w) determined by f{u,v,w) — g{u,v,'w) = h{u,v,w) = 
exists and that this state is stable in absence of diffusion. 

We introduce small perturbations to the steady state as (w,t;,u;) = (u^v^w) -|- ((5m, 5w). Substituting this into 
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Eqs. ([T]), the following linearized differential equations for the perturbations are obtained: 



—Su = fuSu + fySv + fniSw + DuV'^Su, 
at 



dt 



Sv = QuSu + g^Sv + giuSw + Dy'V'^Sv, 



(2) 



—Sw = huSu + hySv + hyjSw + D^^^Sw, 
K at 



where = .^g), fv ~ df ldv\{^u.v,w)^ fw = 9f/9w\(^u,v,w)j ■■■ are partial derivatives at the steady state. The 

following rescaled variables are introduced for convenience: 



Su ~ Su, Sv 



1 k 


Sv, 


Sw = 


k 


1 gu 




hv 



^w, i = V|/„g„|t. 



(3) 



We substitute these variables into Eqs. ([2]) to obtain a set of equations 

( d 

—Sii = moSu + aSv + ?iSw + DuiJ,^^Su, 
dt 

—Sv = j3Su + pqSv + qSw + DyfiV^Sv, 
dt 

—Sw = rSu + ^Sv + sqSw + D^^lV^Sw, 
< dt 



(4) 



where 



mo = 



r = 



fu 



VI/ 








VIA 




K I 


k 




9u 


1 





fw I 

\J\fv^gu\ 



fvdu 



So 



\/\fvgu\ 



\J\fvgu\ 



(5) 



The coefficients a, (3, 7 are determined by signs of /„, gu, h^, a = sign(/t,), /3 = sign(gti), 7 = sign(/it,). The perturba- 
tions {Su,Sv,Sw) are expanded over plane waves as 



(5w = y dkui'') exp -ik-x 
Sv = / dfcy(*^) exp \i^)t-ik ■ X 
Sw = / dkW^^) exp \\^'')i~ik ■ X 



(6) 



where k — {ki,k2,- ■ ■) is the wave vector, and is a growth rate of the plane wave with wave vector k. Thus, we 
obtain the following equations for each wave vector k: 



'mo - Dufik^ 

P 
r 



Po 



1 So 



q 

Duj^J.k'^ 



(7) 



where k = k is the wave number, the magnitude of the wave vector k. Because only the magnitude of the wave 
vector is important, here and below we drop the vector symbols. 
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The condition 

(mo - Duf^k'^ - X^''^ a n \ 

f3 po- D^iik"^ - \(^^ q =0 (8) 

r 7 .so - - A(^V 

should be satisfied for Eq. ^ to have non-trivial solutions. Thus, the linear growth rate A''^-' is determined by the 
characteristic equation 

— {m + p + s)}? + {mp + ps + sm — nr — jq ~ a/3)X — {mps — npr + aqr — ■fmq + f3jn — a(3s) = 0, (9) 

where 

m =m{k) = mo — k^ijDu, 

p=p{k)^Po-k^nDy, (10) 
s =s{k) ~ Sq — k^iiDyj. 

The growth of each plane- wave mode is determined by the real part of A^*^' . The uniform steady state is stable if 
Re(A''^-') is negative for all k. The instability occurs if Re(A('"'^) becomes positive for at least one wave number k = k^. 
Then the uniform steady state is destabilized, leading to spontaneous development of wave patterns with critical wave 
number fcc- If the imaginary part Im( A of the unstable mode is zero, the first critical mode represents a stationary 
plane wave and the Turing instability occurs. On the other hand, if Im(A^'"''^-') ^ 0, the critical mode is oscillatory in 
time and periodic in space, so that the critical mode represents a traveling wave and the wave instability takes place. 

The complex conjugate root theorem tells that a characteristic equation with real coefficients 

\^ +a\^ + h\ + c = Q (11) 

has either three real roots or one real root and a pair of complex conjugate roots. In the latter case, the three roots 
can be written as 

Ai,2=?/'±iw, A3 = (/), (12) 

so that the coefficients are represented as 

a = -(2^- + <?!)), & = + + 2V"/>, c= -{il? +uj'^)(l). (13) 

Combining these three equations, we obtain 

c - a6 = - (lA^ + ^2)0 + (2Vj + </))(V^2 ^^-2 ^ 2V'(/>) 

^2'4^[{i! + 4)f +Lo^]. (14) 

At the threshold of the wave instability, we have = 0, so that c — a6 = 0. At the threshold of the Turing instability, 
we would have = 0, and therefore c = 0. Note that the Turing instability is also possible when the characteristic 
equation has three real roots and one of them becomes positive. It can be easily checked that, also in this case, the 
instability threshold corresponds to c = 0. 

Thus, the wave instability first takes place when, for one wavenumber k ~ kc, the equation 

-^wav ={mps — npr + aqr — jmq + j3^n — a/Ss) + {m + p + s){mp + ps + sm — nr — — a^) = (15) 



becomes satisfied, where m,p and s are given by equations (jlOp with k = k^. 

The Turing instability first takes place when, for one wavenumber k ~ kc, the equation 

1st —rnps — npr + aqr — jmq + /3^n — a/3s ~ (16) 

becomes satisfied, where the coefficients m,p and s are again given by equations ()10p with wavenumber kc. 

It is convenient to introduce the three-dimensional m-p-s space in order to represent these conditions graphically. 
As illustrated in Figures 1 and 2, each of the conditions ([15]) and (fT6)) defines a boundary surface and Eqs. 
determine a straight line F which is parameterized by the wave number k. If the line F touches the boundary surface 
^wav(w,p, s) = or Ist{m,p,s) = 0, then the conditions (|15p or ([T6| are satisfied and the corresponding instability 
takes place. 
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(a) (b) 




-s 



FIG. 1: Boundary surface for the wave instability 7wav(?Ti,p, s) = 0. The hue F touches the surface at a point P. Panel (b) is 
a closeup of (a). Parameters are fixed at n — q = r = a — l3 = ^=l, mo — 0.4, po = 0.4 and so — —1.8. Diffusion constants 
are _D„ = 1, D„ = 1 and D-u, = 14.521. 




FIG. 2: Boundary surface for the Turing instability lBt{m,p, s) — 0. Parameters are fixed atn = g = r = Q = /3 = 7 = l. 



The uniform steady state {u, v, w) should be stable if the diffusion is absent. In terms of of the Jacobian matrix at 
the steady state 

fu fv fu 

J = \ 9u 9v 9w \ , (17) 



hy hit 



this implies that 



det( J) < and tr(J) < 0. 



(18) 
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Thus, the initial point Pq with coordinates (7710,^07*0) on the hne (|10[) should lie inside the stable region in m-p-s 
space. 

The wave instability takes place if the first critical mode is oscillatory. At the threshold of the wave instability, the 
line r defined by Eq. (jlOp should touch the boundary surface /wav = without having intersections with the surface 
1st = 0. 

Suppose that we want to check whether the wave instability can occur when the diffusion constant D^o is varied. 
Let us denote A = Dy/D^ and B = po — moDy/ and consider a plane p = Am + B parallel to the s-axis. The 
line r is always lying on this plane irrespective of D^, because the conditions po = Amo + B and Dy/D^ = A are 
satisfied. A coordinate x is introduced on the plane in such a way that we have x — when m + p = holds. The 
slope of the line F in the x-s space is DuD^/ \/ Dy^ + Dy^ and is nonnegative. The initial point Pq of the line T is 
at (xqjSo) on the plane, where xq = {ttiq +pq)/{A + 1). Because we assume that the necessary conditions (|18|) are 
satisfied, the point (a;o,so) is in the stable region of the x-s space. Two surfaces /^av — and 1st = intersect the 
plane along the boundary curves /wav = and 1st = 0, 



/wav i-^: ^ : A^ B^ /osc 

ist{x,s;A,B) = 1st 



B 



A + r 

B 

A + 1' 



Ax ■ 
Ax 



B 



A + r 

B 



= 0, 
= 0. 



(19) 
(20) 



^+1 

Below we examine the dependences of the boundary curves /wav ~ and 1st = on the parameters 
mo,po, so,n,q,r,a, 13 and 7. This allows us to obtain the parameter conditions under which the wave instability 
can occur. 

Let us examine the shape of the boundary curve /wav = at large values of \s\. If \s\ 1, we can neglect 0{s^) 
terms and obtain 



/wav (a;, s; A, B) PS s [(1 + Afx'^ + s{l + A)x - nr - 79] 
Then, the boundary curve /wav = is given by the equation 



(1 + Afx^ + s(l + A)x -nr --fq = 



(21) 



(22) 



and we have 




{A + iy 



-i± 1 



nr + 7g 



2(^ + 1) 

As follows from Eq. (|2ip . the boundary curve approaches a; = in the limit of large |s| as 



{A + l)xs'^ 



0. 



(23) 



(24) 



Then, the plus sign should be chosen in Eq. ([23| . Thus, the asymptotic boundary curve in the limit s » 1 is the 
hyperbola 



xs 



nr + 7(7 
A + 1 ■ 



(25) 



If the coefficients of this hyperbola is negative, i.e. nr + jq < 0, the boundary curve lies in the region of x > and 
s < 0. In such cases, given that xq > 0, the line F always touches the boundary curve under increasing the diffusion 
constant Dyj, as shown in Figs. [3] and 01 Therefore, we require nr + jq < and > as a part of the sufficient 
conditions. Using definitions of the model parameters ([5]), these requirements can be written in terms of the Jacobian 
matrix as 



fwhu + gu,hy < 0, 
fu + gv > 0. 



(26) 
(27) 



6 




FIG. 3; Boundary curves 7wav = (red) and /at = (blue) when (a) the condition (|30p and (b,c) the condition (|34p are 
satisfied. All species are mobile. Positions of the line F defined by (|10p are shown for two different values of D^]- The boundary 
curve 1st ~ can behave as shown in (b) or (c), depending on the model parameters. 



The boundary curve for the Turing instability Istix-, s; A, B) =0 is given by 

, , _ (A + l)^{Anr + -fq)x -{A + l)^(^7n + aqr) + {A+ l)B{nr - 7g) 
■"^"^^ ^ AiA + iyx"^ -{A-^ -l)Bx-aP{A+iy - B-^ ' ^ ' 

which is a single- valued function of x. 

Let us assume that all three reactants diffuse over the space, Du ^ 0, Dy ^ and ^ 0, which leads to A 0. If 
the denominator on the right hand side of Eq. psp is not zero for all values of cc, then s{x) is a continuous function. 
Figure [3] (a) illustrates the qualitative shape of the boundary curves /wav = and igt = in such a situation. In 
this case, the line F touches the boundary curve /wav = without intersecting /gt = as is increased. As a 
consequence, the wave instability takes place. The conditions are given by 

(A^ - 1)2^2 + ^A{A + 1)2 [a/3(A + 1)2 + < (29) 

which implies 

t f ifuDv + gvDu)'^ ,„„N 

Ugv - fvgu > ^ • (30) 

If the denominator on the right hand side of Eq. ([25]) vanishes at x_ and x+, the boundary curve 1st ~ diverges 
in two ways depending on the values of n, q, r, a, j3 and 7 (Fig. [3] (b) and (c)). In both cases, if xg < a;-, the line F 

{A^-l)B-{A + lf^B^+AapA 
= 2A(A + 1)2 

which is equivalent to 

fu9v - fv9u > (32) 

and 

fuD, + g,Du < (33) 

is satisfied, the wave instability takes place as Dw is increased. 

Next, we assume that the reactant V does not diffuse, that is Dy = 0, leading to A = Q. In this case, the boundary 
curve 1st = is given by 

qx + K{n-qr)+B{nr~q) 
<-)\a=o- Bx + l-B^ ' 
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FIG. 4: Boundary curves Iwav = (red) and 7st = (blue). The reactant V does not diffuse {Dv = 0). The hne F defined 
by (|10p is shown as black lines for two different values of Dw The boundary curve 7st = can behave as shown in (a) or (b) 
depending on the parameters n, q, r, a, /3 and 7. 



which diverges aXx^B— 1/B=pq — 1/po [Figs. 2] (a) and (b)]. In this case, if < — 1/poj the line F touches 
^wav = without intersecting I^t = when is increased. Thus, the instability condition for a system with two 
diffusible reactants is given by 



which implies 



1 

xo<pa , 



gvifudv - fvgu) < 0. 



(35) 



(36) 



On the other hand, if the reactant U does not diffuse, Z?„ = 0, one can choose a plane m = A'p+B' where A' = Du/Dy 
and B' = rriQ — poDu/D^ and derive the condition 



fuifugv - fvgu) < 0. 



(37) 



Thus, sufficient conditions for the wave instability under increasing diffusion constant D^, of the reactant W has 
been derived for the four cases depending on the diffusion constants. The first two conditions ([25]) and ^7} are 
common to all cases. The last condition depends of the diffusion constants, i.e. we have 



fug. - fugu > (^^^ ^ 0, ^ 0), 



fugv - fvgu > and fuDy + g^Du < (Du 7^ 0, Dy ^ 0), 
fuifugv - fvgu) < {Du = 0), 

^ gvifugv - fvgu) < (£)„ = 0). 

These different requirements can however be expressed in a single equation; 

det 



gu 



gv 



+ J^OforanyA<0. 



(38) 

(39) 
(40) 
(41) 

(42) 



with at least two non- vanishing diffusion constants. 

The wave instability may take place also under the variation of other two diffusion constants Dy or Dy. The 
corresponding sufficient conditions for each case can be obtained by permutating three variables u,v and w. 

In summary, the wave instability occurs under increasing diffusion constants Dy or I?.^,, if the conditions 



Qy + tluj > and Qufy + hufw < 



and detf^^'X^^^ , , ^ 7^ for any A < 0, 

\ hy h,„ + D,„A 
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+ fu> and h + fv9u < and det 



fw 




) 



^ for any A < 0, 



(44) 



or 



fu + ,9t) > and fwhu + g^h^ < and det 



( 



fu + DuA 
9u 



fv 

gy + DuA 



^ for any A < 0. 



(45) 



are satisfied. 



We have derived sufficient conditions for the wave instabihty in general three-component reaction-diffusion systems. 
The conditions are formulated in terms of the Jacobian matrix elements at a steady state and of the diffusion constants 
They do not depend on model details. Once these conditions arc satisfied, the wave instability occurs as we increase 
the diffusion mobility of one of the reacting species. 

Our general results are applicable for systems of various origins, including biological, chemical, physical and ecolog- 
ical systems. Our analysis has revealed that the wave instability may occur even if one of three rcactants is immobile. 
This result can be important in a variety of applications involving both diffusible and non-diffusible reactants. 
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